WIMPs in a 3-3-1 model with heavy sterile neutrinos 
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Abstract 

In this work we build a gauge model based on the SU(3) C ® SU(3)l ® U(1)n symmetry with heavy 
neutrinos and show that we can have two weakly interacting cold dark matter candidates in its spectrum. 
This is achieved by noticing that a global U(l) symmetry can be imposed on the model in such a way that 
the stability of the dark matter is guaranteed. We obtain their relic abundance and analyze their compatibility 
with recent direct detection experiments, also exploring the possibility of explaining the two events reported 
by CDMS-II. An interesting outcome of this 3-3-1 model, concerning direct detection of these WIMPs, is a 
strong bound on the symmetry breaking scale, which imposes it to be above 3 TeV. 
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I. INTRODUCTION 



The Dark Matter (DM) problem constitutes a key problem at the interface among particle 
physics, astrophysics and cosmology. The observational data accumulated in the last decade point 
to the existence of a non-negotiable amount of non-baryonic DM, whose identity is still unknown. 
Since the standard model of electroweak interactions (SM) does not provide any candidate for such 
invisible component of matter, this problem is an indication for physics beyond the SM. Different 
measurements coming from cosmic microwave background radiation (CMBR) [HJ], galaxy rota- 
tion curves 0], gravitational lensing [|[] and structure formation 0| etc., confirm that, besides its 
undoubted evidence, it must contribute to around 22% of the total energy density of the Universe. 
Nowadays it is known, due to numerical simulations to reproduce the structure formation [|5|], that 
the matter present in our Universe is dominated by cold dark matter (CDM), and precise measure- 
ments of its relic abundance impose strong constraints on various new physics models. However, 
the relic abundance alone is not enough to point out all the properties of the CDM, even already 
assuming that the CDM is represented by a weakly interacting massive particle (WIMP). Col- 
liders, and particularly the Large Hadron Collider (LHC), have the potential for discovering and 
identifying new particles not predicted by the SM. By the other side, most of the information about 
the nature of the CDM might be extracted from existing direct detection experiments, through the 
scattering of WIMPs with nuclei [|5p and therefore we will focus on them. 

The WIMPs are the most studied CDM candidates and arise naturally in several theoretical 
frameworks such as Supersymmetry f7|, §], Universal Extra Dimensions [@], Little Higgs Mod- 
els [[10p , Technicolor [fTT|], etc., but since all these theories remain hypothetical [JT^], it is equally 
worthwhile to tackle less conventional possibilities. For this reason, we are going to explore a 
small gauge extension of the electroweak sector of SM, SU(3)c ® SU(3) L ® U(1) N , 3-3-1 for 
short (for a nice review see Ref. [|T3l]). This extension can be accomplished by a class of mod- 
els [ |14[ ] that have intriguing features such as: the models are anomaly free only if the number of 
families is a multiple of three allied to the condition of QCD asymptotic freedom [|15J]; the electric 
charge quantization and the explanation of the vector-like nature of the electromagnetism are natu- 
rally achieved in the absence of anomalies [JT^]; there is room for lepton number violation JT7| ] and 



new sources of CP violation Ql8Q , crucial features to approach baryogenesis and/or leptogenesis, 
among other nice characteristics of the model. 

The CDM problem in 3-3-1 models was already studied in different situations, with self- 



interacting DM |]T9[], a scalar bilepton (a particle that carries two units of lepton number) WIMP pQ[] 



and a supersymmetric self-interacting DM []2T|]. Here we are going to consider a variation of the 
model of Ref. p0| ] and make a deep analysis of the CDM candidates and their respective abun- 
dance and direct detection. In the work of Ref. [ |2D[ ] there are light right-handed neutrinos in the 
triplet representation of SU (3) L and no singlet neutrino, while the version of the model developed 
here contains new left-handed neutrinos in the fundamental representation of SU(3) L instead of 
right-handed neutrinos. For this reason we will call it 3-3-1LHN for short. This 3-3-1LHN model 
was inspired by the first attempts to enlarge the electroweak gauge symmetry SU L {2) ® U Y (1) 
to SU L (3) ® U N (1) p2|]. It is amazing that a global £7(1) symmetry can be imposed that not 
only simplifies the Yukawa Lagrangian and the scalar potential but also stabilizes the lightest of 
the new particles charged under this symmetry, providing candidates for explaining the CDM 
problem. Our main goal is to get WIMP CDM candidates from the spectrum of 3-3-1LHN in 
agreement with most recent bounds from direct detection experiments, namely CDMS-II [ ]23[ ] and 



XENON [[24J, and investigate the region of parameter space which is well suited for explaining 



the positive signals observed by the CDMS-II experiment Q25Q. We do not take into account 
DAMA [ |26[ ] and CoGeNT [ p7| ] results in this work since our model does not have any allowed 
region in the parameter space that explains the observed signals by these experiments. Neverthe- 
less, we remark that the parameter space favored by these experiments is mostly in conflict with 
all other detection experiments, although it remains an intriguing challenge to be solved. 

We start by briefly describing the model in section ||, introducing its main ingredients. In 
section |l| we discuss the relic abundance computation as well as the direct detection bounds for 
the CDM candidates and analyze the compatibility of our model with the positive signal from 
CDMS-II. We present our conclusions in section [IV|. 

II. THE 3-3-1LHN MODEL 

In the 3-3-1LHN model the leptons are accommodated in triplet and singlet representations as 
follows (we indicate the SU (3)c ® SU (3)l ® U (1) n transformation properties in parentheses), 



faL 



( v \ 

a, 



\N a J 



(1,3, -1/3) , e aR ~ (1, 1, -1) , N aR ~ (1, 1, 0), (1) 



where a = 1, 2, 3 represents the family index for the usual three generation of leptons, while 
N a (L,R) are new heavy neutrinos representing new degrees of freedom in this model, and it is 
this assumption that makes the 3-3-1LHN model substantially different from the proposal studied 
before 0]. 

In the Hadronic sector, the first generation comes in the triplet representation and the other two 
are in an anti-triplet representation of SU L (3), as a requirement for anomaly cancellation. They 
are given by, 
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(2) 



where the index i = 1, 2 where chosen to represent the first two generations. The primed quarks 
are new heavy quarks with the usual fractional electric charges. 

In order to generate SM fermion masses, three scalar triplets are introduced, 
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with Tj and x both transforming as (1 , 3 , —1/3) and p transforming as (1 , 3 , 2/3). 
In general, a discrete Z 2 symmetry is usually assumed, transforming the fields as, 

(X, P, e aR, N aR , u aR , d' iR , Q 3L ) - (x , p , e aR , N aR , u aR , d' iR , Q 3L ) 

which leads to an economical model with a simplified Yukawa Lagrangian Q, 

- C Y = fijQiLX*d' jR + hzQzLXU 3R + giaQiLV*daR 
+ h 3 aQ3LV u aR + 93aQ3Lpd a R + h ia Q iL p*U aR 

+G ab f aL pe bR + g' a JaLX N bR + H.c. 



(4) 



(5) 



1 We are going to see that this discrete symmetry can be replaced by a global U(l) symmetry, with the effect of 
producing the same terms in the Lagrangian of the model, but also advantageous in stabilizing our CDM candidates. 
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Again, in these expressions we are using the family indices i = 1, 2 and a = 1,2,3. 

The most general scalar potential that we can construct which obeys the above discrete sym- 
metry has the form, 

v(v, P, x) = p 2 x x 2 + p 2 v v 2 + p 2 p p 2 + AiX 4 + W + A 3 p 4 + 

A 4 (x t x)(^ t ^) + AsO^xXpV) + h{rfv)(p*P) + 

ArCx^X^x) + A 8 (x t p)(p t x) + A 9 (77 t p)(P t ^) 

-^€«ViXib + H.c. (6) 

It is well known that this potential is appropriate to induce the desired spontaneous symmetry 
breaking pattern of the electro weak gauge symmetry, SU(3)l ® U(1)n to SU(2) L ® U(l)y and 
finally to U(1)qed 9 generating the masses of gauge bosons and fermions through the Higgs mech- 
anism. 

We also write the currents involving the non-hermitian vector bosons for leptons and quarks, 
since it is going to be necessary in observing the existence of an extra global symmetry in the 
3-3-1LHN model. It reads (see the fourth paper in Ref. [|14|]), 

C NH = -X [i^felW; + Nl*fe a L V+ + v a L ^N a L U° 

+ (usL-fu'u - d[ L ^d iL ) Ul + h.c] , (7) 

where we have defined W+ = ~^(W^ - iW%), as usual, V~ = ^(Wj - iWfi and C/J = 
75 (W 7 ^ ~~ iW^). The three remaining neutral gauge bosons, A^, and Z'^ couple to the fermions 
in a diagonal basis and do not influence the discussion on the new symmetry that follows, so we 
do not present their currents here. 

Now we notice that there exists a new extra global symmetry in this model, which we call 
U (1)g, with the following assignments of G charges carried exclusively by the 3-3-1 model fields, 

G(N L/R , u' 3L/R , d' iL/R , V-, Ul x°, X", v'°\ P") = +1 • (8) 

All the other fields transform trivially under this symmetry. At this point we notice that we could 
have started with this U(1)g global symmetry from the beginning, without imposing the previously 
mentioned discrete symmetry, Eq. (||). In other words, if we replace the Z 2 global symmetry by 
the U(l )q global symmetry, we recover the same Lagrangian terms as given in Eq. @ and Eq. @ 
with no new term to be added. The advantage of this continuous symmetry is that the G charged 



fields (we call them G-fields for short) always appear in pairs, guaranteeing that the lightest one 
is stable. Next we identify the mass eigenstates in the 3-3-1LHN model such as to select which 
neutral G-fields can be a potential CDM candidate []. 



A. The mass eigenstates 

In order to achieve spontaneous symmetry breaking, we suppose that the neutral scalars 
(jf-> P°> develop a vacuum expectation value (VEV) according to, 

V°, P°, x'° -> ^K P ,x' + Rn,P,x' + U v,P,x f ) , ( 9 ) 
where we make the reasonable and simplifying assumption that the remaining neutral scalars 
(?/ , x°) do not develop VEVs g. 

From this pattern of symmetry breaking, we observe that the U(1)g symmetry forbids Majorana 
mass terms for the neutrinos and no mixing appears among the new neutrinos with the standard 
ones. This turns them into truly sterile Dirac neutrinos. Moreover, for sake of simplicity, we 
consider that the mass matrix of the charged leptons, new neutrinos and of the new quarks all 
come in diagonal mass bases with normal hierarchy. 

Considering the vacuum structure in Eq. @, the mass matrix of the new neutrinos and quarks 
take the form, 

M Na = ^v x , , (10) 

and 

L 

respectively. If we assume that g' u < g' 22 < g'33 in the first of these equations, the lightest heavy 
neutrino is identified with Ni. Regarding the standard neutrinos, we assume here that their tiny 
masses are due to effective dimension-5 operators as first implemented in Ref. 

As for the scalar mass matrices, we first need the minimum conditions from the potential in 
Eq. @, given by, 



^ = i i«v> (id 



2 We present the trilinear couplings for the G-fields in appendix A. 

3 If we take non-trivial VEVs for these scalars we would still obtain the complete mass spectrum of the model with 
only additional complexity in the mixing of gauge bosons and scalars. However, this would also break the U(1)g 
global symmetry, yielding an unwanted Goldstone boson in the spectrum. 
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^ + A 3 ^ + f4 + ^-^ = 0. (12) 

Although the trilinear coupling / in Eq. @ is a free mass parameter, in this work we make 
the assumption that / is of the order of the 3-3-1 symmetry breaking scale, v x * 9 supposed to be at 
TeV scale, while v p and v v (<C v x t) have to be at the electroweak breaking scale, v « 246 GeV, 
since they fix the Z and W ± gauge boson masses [fT4|], being related by v*-\-v 2 = v 2 . We then 
choose f — V ~y and v p = v v = just to simplify the diagonalization procedure of the scalar 
mass matrices. 

Substituting the Eqs. (Q) and ( |T2| ) into the scalar potential, Eq. @, we can obtain the mass 
matrices for the neutral scalars in three different bases, a scalar, (R x > , R v , R p ) 9 a pseudo-scalar 
one, (I x f , , I p ) 9 and a complex scalar basis, (x ot , ^ /0 ). 

Since no fine-tuning is assumed we can take some simplifying relations here in order to obtain 
the mass eigenstates. Namely, A 4 = A 5 = 0.25 and A 2 = A 3 . This assumption, which in principle 
would demand some kind of symmetry to guarantee such equalities, may have some implication 
on the mixing of interaction eigenstates which could change somehow our results. However, at 
this point we still do not have a consistent way of considering more general scenarios where this 
diagonalization can be numerically implemented, an issue we hope to develop in the future. For 
this reason we are going to perform our computations in this framework, keeping in mind that 
some different outcome could emerge in a more general scheme, which would have our scenario 
as a subset of possibilities. 

We then find the following mass eigenvectors in the basis (R x > , R v , R p ), 



x j 



S 1 = R 



with respective mass eigenvalues Q, 



2 



4 We notice that the real and pseudo-scalar mass eigenvalues in Ref . pQ| ] are lacking a factor of two, while the WIMP 
complex scalar has a correct factor. This does not change the qualitative results and conclusions in that work, 
although tiny quantitative corrections are implied wherever the Higgs boson play some role. Here we took those 
missing factors into account. 
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M 2 H =v 2 (2X 2 + X 6 ). (14) 



In the pseudo-scalar basis, (I x > , I v , I p ) we find the mass eigenstates, 
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where, 7f and l£, correspond to Goldstone bosons and P\ is a massive pseudo-scalar that remains 
in the spectrum whose mass is, 

< = ^ + y). (16) 
Also, in the basis of complex neutral scalars, (x° , we g et the mass eigenstates, 

r — Vx ' a/ _l _ n /0 * 



*v«/l + £ Jl + 5 2 



X' 



where is recognized as the Goldstone boson eaten by the gauge bosons U° and U°* and 4> has 
a mass, 



Mj = ^% 2 + ^]. (18) 

Considering the two bases of charged scalars, (x~ •> P /_ ) an d {rf~ , p~), we obtain the following 
mass eigenstates, 



= ^fa~ + P~), (19) 



which can be checked to be the same eigenvectors as in Ref. QZOfl when we take the limit v x > 3> v. 
Their mass eigenvalues are, 

M 2 = + A 9 v 2 . (20) 
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The remaining eigenvectors are two Goldstone bosons given by, 



h 3 = 1 -- (X ~ -JV ) , 



K = ^(v~-p-)- (2D 

Finally, from the gauge invariant scalar kinetic terms (not shown here) and using Eq. (|^), we 
easily obtain the gauge boson masses [jPJ], 

2 1 2 2 



m z = m w±/%» 
1 

4(3-4? w ) L ~ w "x' ' c 2 w ' c? w 



2 2 1 2/ 2 , 2\ 

m v ± = m u0 = -g (v x , + v ) , 

m z> = Ti~o~Tj2~\ l^wvl + o - ] , (22) 



where we have defined the Weinberg mixing angle through sin 9 W = s w (as well as cos 9 W = c w ). 
Notice that we have neglected the mixing between the neutral gauge bosons Z and Z' , which is 
constrained to be very small (see the fourth paper in Ref. [|14|]). 

With all the mass eigenstates identified as above we are able to consider the stability of the 
neutral G-fields. In the 3-3-1 model with right-handed neutrinos studied in Ref. PQ], cj) was the 
same combination of interacting neutral scalar G-fields as in this 3-3-1LHN model, but there this 
scalar carried two units of lepton number instead. There is no other neutral G-field scalar in the 3- 
3-1LHN model, and the only neutral G-field vector boson is U°. These, together with the lightest 
heavy G-field neutrino, Ni, are the potential CDM candidates of this model, although they cannot 
be simultaneous candidates since they couple to each other plus some standard model particle, as 
explicitly shown for the trilinear couplings in appendix A. Thus, it is enough to make one of them 
the lightest particle among the G-fields, which provides a stable CDM candidate. We would like 
to stress that although the gauge boson, [7°, could be the stable G-field, it leads to a too much 
suppressed relic abundance, thus we do not consider it henceforth as a third CDM candidate. 

Having found the mass spectrum of the 3-3-1LHN model we identify cp and Ni as our possible 
DM candidates (again, U° is a candidate too, but extremely underabundant) by enforcing that 
one of them be the lightest G-field, we will next determine their relic abundance and analyze the 
WMAP favored parameter space region under direct detection experiments. 



III. RELIC ABUNDANCE AND DIRECT DETECTION 



Among the CDM candidates, the WIMPs are the most intriguing ones since their thermal cross 
section, which is roughly at the electroweak scale, naturally leads to the appropriate relic density. 
The scenario goes as follows: a WIMP which is in thermodynamic equilibrium with the plasma 
in the early Universe decouples when its interaction rate drops below the expansion rate of the 
Universe. In this way we have first to check that the CDM candidate besides being stable (or 
meta- stable), either freezes with the right relic abundance [p]] or, at least, represents the majority 
of CDM constituting a subdominant scenario. Secondly, since nowadays we have some direct 
detection experiments available [^3|, ^4p, it would be desirable that our candidate has at least some 
chance of being detected in the near future or, more remarkably, to explain positive signals such 
as the events in excess observed by CDMSII. 

First we will describe the computational procedure used to get the relic abundance of the 3-3- 
1LHN CDM candidates, cp and Nu and present some scatter plots showing our results for different 
regimes of the parameter space. Lastly, we will discuss a little bit about the direct detection method 
and compute the WIMP-nucleon cross section of our candidates and investigate its feasibility in 
light of CDMS and XENON bounds and also the possibility of explaining the recent CDMSII 
signal. 



A. Relic Abundance 

In order to obtain the WIMP abundance in its decoupling stage we need to solve the Boltzmann 
equation which gives the evolution of the abundance of a generic species in the Universe as a 
function of the temperature, 

where, g* is the effective number of degrees of freedom available at the freeze-out temperature, 
M p is the Planck mass, Y is the thermal abundance or number density over entropy (while Y eq is 
the abundance at the equilibrium epoch) and (av) is the thermal averaged cross section for WIMP 
annihilation times the relative velocity. The particle physics information of the model enters in 
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this cross-section which includes all annihilation and co-annihilation channels, 



J2m f N9 ^^(-Ip^E^;^ 8 ) 

GV ) = ^7v^ 2^7 7^ > (24) 

2T(^g i m i K 2 {m i /T)) 



where ^ is the number of degrees of freedom that characterizes the species involved, a^-hi the total 
cross-section for annihilation of a pair of particles with masses m*, rrij into some SM particles 
(fc, /) with respective masses and m/, while p^- is the momentum of incoming particles in their 
center of momentum frame. 

The relic density is obtained by integrating from T = oc to T = T where T is the temperature 
of the Universe today, yielding, 

nh 2 = 2.742 x 10* MwiMP Y(T ) (25) 
GeV v J 



Our results are obtained by using the package micrOMEGAs [|29|], which computes this relic den- 
sity numerically for a given model. The task would reveal unfeasible analytically since many 
interactions participate in the annihilation process at freeze-out. We have also implemented the 
3-3-1LHN model in the package lanHEP [ ]30| ] that furnishes the model files to be used in mi- 
crOMEGAs, making the task of computing the relic density much easier and reliable. The most 
significant processes which contribute to the abundance of our CDM candidates, Ni and 0, sepa- 
rately, are shown in Figs. [I]-^. After using the procedure described above we then show the results 



Pi 



Z' ^ Z' Z' ^v + ,u° 

>wv< >wv< >vw^ 

FIG. 1. The main processes which contribute to the abundance of N\, where I = e^[x^T^n e ^n yb ^n T and 
q — u^d^ c, 5, 6. 

for each candidate. In Fig. |3] we show the relic abundance for the heavy neutrino WIMP, Ni, for 
v x i = 3 TeV and 4 TeV. We should remark that the masses of Z' and Pi depend only on the values 
of the VEVs, and will not change as we vary the several coupling constants in the model, while 
the Si mass contains an additional free coupling constant, Ai. We then vary the Si mass instead 
of Ai, and the range considered in this case is 400 GeV < M Sl < 4.5 TeV for v x t = 3 TeV and 
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FIG. 2. The main processes which contribute to the abundance of (j) 

600 GeV < M Sl < 6 TeV for v x > = 4 TeV. Hence, the only relevant varying parameter (besides 
the neutrino mass) is M Sl , which leads to a denser region in the abundance for large sterile neu- 
trino masses. The region in accordance with WMAP7, 0.098 < Qh 2 < 0.122, is shown between 
the red bars. One can see that a change on v x > is not going to affect appreciably the shape of the 
abundance, while it considerably change its quantitative aspect, diminishing the favored WMAP7 
region for lower values of v x > . 
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FIG. 3. Relic abundance for the heavy neutrino N\ with the region in accordance with WMAP7, 0.098 < 
tth 2 < 0.122, shown between the red bars. We used 400 GeV < M S i < 4.5 TeV and v x > = 3 TeV in the 
left panel and 600 GeV < Msi < 6 TeV and v x > = 4 TeV in the right one. 

As for the scalar 0, using the same arguments used earlier for Nu we observe that the only 
parameters which control its abundance are the cj) mass and the masses of the Higgs and the scalar 
S2. Let us remark that the S 2 mass depends on the same couplings as the Higgs mass and can be 
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considered to be constant, since in the range of Higgs mass employed in this work, 115 GeV to 
300 GeV, the S 2 mass change only about 5 GeV. Hence, the abundance of cj) is generally governed 
only by the Higgs mass. Nonetheless, the first process in Fig. ^ may be the most relevant when the 
produced quarks are heavy, and then we also vary the intermediate exotic quark mass parameter 
in the range 600 GeV< M q t < 2 TeV In order to see the effect of varying the Higgs mass we 
show two plots in Fig. |] containing our results for the abundance with M H = 115 GeV and 
M H = 300 GeV. Comparing the two panels we conclude that the abundance of (ft is considerably 




200 400 600 800 1000 200 400 600 800 1000 



M (GeV) M* (GeV) 

FIG. 4. The abundance of the scalar (f) for two distinct values of the Higgs mass. The left panel is the 
abundance for M H = 115 GeV and the right one is for M H = 300 GeV. We used v x > = 3 TeV. 

modified by the Higgs mass, with a light Higgs boson offering a denser region on the parameter 
space. 

In summary, the model contains two interesting CDM candidates in two distinct regimes: one 
where Nu a sterile neutrino, can account for the whole CDM and another where 0, a scalar, is the 
CDM. Both can be stable (but not simultaneously, unless they are degenerate) thanks to a global 
U(1)g symmetry, under which only some of the new particles are charged, implying that they 
are always produced in pairs, which resembles something like an R-parity, though it is related 
to the a continuous symmetry instead. Once observed that our candidates can account for the 
total CDM abundance, we need to check if they are in agreement with the last constraints from 
direct detection experiments and through this condition, we are going to assess a constraint on the 
symmetry breaking scale of the model. 

We also want to check if there is some room to explain some of recent claims of a light CDM 
positive signal in CDMSII, which may be possible for the scalar 0, whose mass can be made 
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naturally small. 



B. Direct Detection 

After their decoupling, the WIMPs can cluster and form a local density of CDM surrounding 
us. Therefore the space at Earth location is supposed to be permeated by a flux of these particles 
characterized by a density and velocity distribution that depend on the details of the galactic halo 
model. If these WIMPs are allowed to interact with nuclei, through more fundamental interactions 
with quarks (for a good review see []32p), then it is possible to directly detect them by measuring 
the recoil energy (Q) deposited in the detector material, given by, 

Q = 2^-, (26) 

m N 

where \i r M w m N j \M W + m N ) is the reduced WIMP-nucleus mass, M w the WIMP mass, m N 
the nucleus mass and v is the minimal incoming velocity of an incident WIMP. 

Measuring the energy deposited by the WIMP and making some assumptions about the halo 
model, we can infer the spin independent WIMP-nucleus cross section at zero momentum transfer, 



using the standard procedure described in [g, [29|, [32|] 



a = M ( Z f p + (A - Z)f n f . (27) 

7T 

where Z is the atomic number, A is the atomic mass and f p and f n are effective couplings with 
protons and neutrons, respectively, and depends of the particle physics input of a given model. 



Since the DM experiments such as the Cryogenic DM Search (CDMS) []23|] and the liquid noble 



gas XENON P4 | contain nuclei with different atomic masses, its useful to define what we call the 



WIMP-nucleon cross section when f p = f n , 



,2 



where is the WIMP-proton/neutron reduced mass. The assumption f p = f n is valid for most 
models, but there will be instances in our model where this fails to be true, as we will point later 
for the case of Ni. 

These experiments have been trying to observe WIMP events, but in most of the cases no event 
have been detected and hence they were able to impose strong limits in the WIMP-nucleon cross 
section instead. Nevertheless, recently the CDMS collaboration has reported its results of the final 
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data runs of the CDMSII and observed that two candidate events have survived after application 
of many discrimination procedures. The probability of observing two or more background events 
is 23%, which means that the two events neither provide a statistically significant evidence for 
CDM, nor can be rejected as background. Many works have been done interpreting these two 
candidate events as WIMP signals in different frameworks p3p. Here we will first investigate if 
the 3-3-1LHN CDM candidates satisfy the bounds imposed by these experiments and also search 
for an explanation to the events in excess observed by CDMSII. 

The process which contribute to spin independent cross section of Ni and cj) are shown below, 
in Fig. H and |^. Well, after discussing a little bit about the direct detection method and showing the 
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FIG. 5. Processes which contribute to the WIMP-nucleon cross section of N\ 
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FIG. 6. Processes which contribute to the WIMP-nucleon cross section of (j). 

processes which contribute to the WIMP-nucleon cross section we are able to show and analyze 
the results for each candidate. 

The scattering processes of Ni with quarks are exhibited in the Fig. [5|. Using the fact that the 
vertices involving the gauge boson Z' possess only gauge couplings and that the scalar Pi couples 
to Ni proportionally to its mass, the only free parameters related to the WIMP-nucleon cross 
section of Ni are its own mass and v x >. To see how our results are modified by the value of iy we 
evaluate the spin independent WIMP-nucleon cross section at zero momentum transfer limit given 
in Eq. (@), for iy varying from 2 TeV to 4 TeV, which we present in the Fig. |7|. Actually, in the 
case of TVi, the WIMP-nucleon coupling with protons is about one order of magnitude higher than 
the coupling with neutrons, and we choose to plot the WIMP-proton cross section since it is more 
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strongly constrained than the neutron one in this case. 
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FIG. 7. The WIMP-proton cross section for N\. From top to bottom, the curves represent the variation of 



v x t in the range 2 TeV< iy < 4 TeV. The data used in the exclusion curves were obtained using [34]. 



From the Fig. [7] we might realize that the heavy neutrino constitutes a nice CDM candidate 
obeying the most recent bounds from direct detection experiments if iy > 3 TeV. The changing 
in the scale of symmetry breaking shows us that raising the values of v x t we make the model 
safer if the experiment sensitivity grows. This is an interesting result for this model because the 
direct detection experiments are strongly constraining the breakdown of 3-3-1LHN symmetry to 
be above 3 TeV. The gap in the results on this figure appears because it refers to the overabundant 
regime (Qh 2 > 0.122) whose points were not included in the plot. 

Finally for the scalar 0, we can also calculate the WIMP-nucleon cross section taking into 
account the possible processes shown in the Fig. To understand how many free parameters 
are really important to the WIMP-nucleon cross section we provide some details in what follows. 
Since the exotic quark Yukawa couplings and scalar couplings are naturally of the order one, the 
cross section dependence on them can be translated into their masses, while for the S 2 scalar this 
reflects directly on v x t, since its mass is given by Eq. ([141). Also, the gauge boson Z' contribution 
involves only gauge couplings. Therefore, the only free parameters are the exotic quarks, Higgs 
and cj) masses, besides v x >. To precise how much the Higgs mass is important to our results we 
show the WIMP-nucleon cross section of (j) for different Higgs masses in Fig. |[ 

Comparing the first three plots in Fig. |8] we conclude that the WIMP-nucleon cross section of 
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FIG. 8. The WIMP-nucleon cross section for (j). The first panel is for Mh = 115 GeV. The second one is 
for M H = 156 GeV. The third panel is for M H = 300 GeV. We used v^> = 3 TeV in the first three plots 
while this parameter is varying and Mh = 156 GeV in the fourth plot, where the darker colors indicate 
those points in agreement with WMAP7, while lighter colors represent the region below the WMAP7 upper 
bound. The data used in the exclusion curves were obtained using [|34|]. 



the scalar (ft is very sensitive to changes in the Higgs mass and that the best parameter space in 
agreement with WMAP constraints is obtained for a light Higgs boson, in particular for a mass 
around 150 GeV. The spread points in each plot beyond w 500 GeV are due to the changing 
in the masses of the exotic quarks in the range 636 GeV< M q / < 2 TeV, which does not affect 
the cross section for lower WIMP masses. The cross section dependence on the S 2 mass, which 
is basically v x t 9 have an impact on the results as shown in the fourth panel of Fig. |8|, where we 
exhibit the WIMP-nucleon cross section behavior for different values of v x » and M H = 156 GeV. 
In that plot the points in lighter colors represent a region of the parameter space corresponding to 
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fi/i 2 < 0.122, while the darker colors are the region in agreement with WMAP. We do this with 



the purpose of showing that direct detection bounds on the scalar (f) seems to disfavor v x > below 
3 TeV, as in the case of the sterile neutrino N\. It is noticeable that most points for v x > = 3 TeV 
are ruled out for large WIMP masses, and a light (j) is favored, while for v x > = 4 TeV the model 
is on the verge to be tested for the whole (j) mass range. This reasonably high symmetry breaking 
scale is something to be taken into account when looking for 3-3-1LHN signals at LHC, a task we 
intend to perform soon. 

In brief, we have checked that the two CDM candidates separately satisfy the exclusion limits 
from the most restrictive DM detection experiments. Now we will show that the model also pro- 



vides an explanation to the two excess events observed by CDMSII Collaboration Q23Q. Since the 
scalar cp is the only candidate which can have low mass in agreement with these limits, it is the 
only one capable of representing those excess events. Computing again the WIMP-nucleon cross 
section only for low masses letting the Higgs mass free to vary from 115 GeV to 300 GeV, we 
obtain the behavior depicted in Fig. [9|, We can observe that for the scalar 0, there exists a region of 
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FIG. 9. The WIMP-nucleon cross section for low masses of (f) with 115 GeV < Mh < 300 GeV. We used 



3 TeV. The data used in the exclusion curves were obtained using [34]. 



parameter space with M$ below 60 GeV for which 10 43 cm 2 < a si < 10 44 cm 2 , that reproduces 



the two candidate events reported by CDMSII [[23J] and is not excluded by the recent bounds from 
XENON 100. In addition to this, comparing the parameter space for distinct values of the Higgs 
mass we conclude that if these events are really a WIMP signal due to the scalar 0, it prefers a 
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Higgs boson with M H ^150 GeV (see Fig. |8]) and in this way the solution of the DM problem 
brings to us some hints on Higgs search. It is important to say that the results shown above were 
obtained with only two free parameters, which are the scalar potential coupling constant, A 7 , and 
the Higgs mass. 

Finally, we should mention that this work has not only enlarged the possibilities of candidates in 
the 3-3-1 model, as compared to Ref. [ l20[] , but also a deeper analysis was carried out considering a 
wider range of parameter space. This was achieved by the introduction of a heavy neutrino into the 
spectrum allowing for a new global U(1)g symmetry that would no be possible in the 3-3-1 model 
with right-handed neutrinos [J20|]. This new symmetry was crucial to establish that the lightest new 
particles charged under U(1)g are stable. 

IV. CONCLUSIONS 

We have studied a 3-3-1 model with heavy sterile neutrinos and observed that the model ac- 
commodates a new extra global U(1)g symmetry, that makes possible the identification of three 
CDM candidates in its mass spectrum. One of them, a non-hermitian vector boson U° was not 
considered in our analysis because it does not provide enough amount of CDM. The remaining 
two are a neutral scalar, 0, studied before in another version of the 3-3-1 model []20| ] in a very 
restricted scenario (and stable thanks to a lepton number symmetry), and the lightest of the heavy 
neutrinos, Ni. We have shown that the scalar (j) and the sterile neutrino N\ can account for the 
total CDM in agreement with WMAP7 data. We then computed the scattering cross section of 
our WIMPs with nucleons, in order to comply with recent direct detection experiments, CDMSII 
and XENON 100, and concluded that there is a large range of the parameter space that obeys their 
exclusion limits. For the scenario where (f) is the WIMP DM we also concluded that a Higgs mass 
of about 150 GeV is favored by these limits. 

Besides, an interesting outcome has emerged from our analysis concerning the direct detection 
of our CDM WIMPs, Ni and 0, which is the fact the characteristic symmetry breaking scale of 
the 3-3-1RHN model, v x * 9 has to be larger than about 3 TeV so as to evade the current exclusion 
limits from CDMSII and XENON 100. This is an important feature to be considered in testing this 
model at LHC (a work to be developed elsewhere). 

Finally, we saw that the scalar cj) might reproduce the two excess events reported by CDMSII in 
the Fig. |9|. Our results imply that the model can either satisfy the exclusion limits and/or explain 
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the positive signal observed by CDMSII, pointing to a Higgs mass below 300 GeV with strong 
bounds on the 3-3-1RHN symmetry breakdown scale. 
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APPENDIX A 



In the tables below we present all the triple interactions and couplings involving the G-fields in 
the 3-3-1LHN model, relevant to determine the stability of our CDM candidates. Here we define, 
e as the electric charge and, 

Ri = ^1 + ^, #2 = ^2 + ^-, g w = l- 2s 2 w , a, = 3 - 4s 2 w 

For simplicity, no hermitian conjugate interaction is included in these tables and interactions 
that are already included in one table are not present in the others. 
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